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Our point of departure are the unitary dynamics of closed quantum systems as generated from 
the Schrodinger equation. We focus on a class of quantum models that typically exhibit roughly ex- 
ponential relaxation of some observable within this framework. Furthermore, we focus on pure state 
evolutions. An entropy in accord with Jaynes principle is defined on the basis of the quantum ex- 
pectation value of the above observable. It is demonstrated that the resulting deterministic entropy 
dynamics are in a sense in accord with a transient fluctuation theorem. Moreover, we demonstrate 
that the dynamics of the expectation value are describable in terms of an Ornstein-Uhlenbeck pro- 
cess. These findings are demonstrated numerically and supported by analytical considerations based 
on quantum typicality. 
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1. Introduction Fluctuation theorems as general rules 
controlling the entropy production in all sorts of phys- 
ical systems have been in the focus of nonequilibrium 
physics for roughly two decades [iHl]- Especially fluc- 
tuation theorems describing distributions of work have 
been extensively analyzed for classical deterministic and 
stochastic dynamics [5|-|7[. In the context of quantum 
mechanics most investigations in this field are also on 
quantum work fluctuation relations (st-fllj. These ap- 
proaches usually require some notion of work performed 
by a time-dependent Hamiltonian and the system to be 
initially in a canonical Gibbs state [l2| . In this Letter in 
contrary we do not consider any work but fluctuations of 
entropy in a class of closed non-driven quantum systems. 
Furthermore, we focus on pure states rather than Gibbs 
states. The entropy will be given below simply as a func- 
tion of the (pure) state of the system. Thus "entropy 
fluctuations" here refer to deterministic but irregular ap- 
pearing (small) deviations of the entropy evolution from 
a smooth mean behavior. To some extent comparable 
scenarios have been studied in However, none 

of these works considers the specific entropy definition 
used in this paper in combination with the time evolu- 
tion of pure states directly arising from the Schrodinger 
equation. 

In general the fluctuation theorem (FT) is said to hold 
if the probability, or in the case at hand rather the rel- 
ative frequency, of mean entropy productions (averaged 
over a time step r) E^., obeys the following relation: 



(1) 



We analyze relaxation processes in non-driven systems, 
thus we do not allude to time reversed trajectories as 
often done in this context, i.e., here ([T|) describes the 
probabilities of producing/annihilating certain amounts 
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of entropy along a time-evolving trajectory. Moreover 
Eq. ([T]) is meant as a transient FT ([T|) is meant to apply 
for time steps r that are short compared to the time scale 
of the relaxation dynamics (t <^ tj^). 

The probability density function does not necessarily 
have to be Gaussian to fulflll the FT. But if it is Gaus- 
sian which is (approximately) the case in our numerical 
calculations (see below), it has to take on the following 
form in order to fulflll the FT: 
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with Eq being the mean entropy production and the 
variance. 

2. Introduction of the Model and Expectation Value 
Dynamics Our numerics are based on a model which is 
designed to represent the most simple closed quantum 
model featuring exponential relaxation of an expectation 
value (a similar model and its overall behavior have al- 
ready been established in [13 )• The model, which is 
sketched in Fig. [1] describes a "quantum two site hop- 
ping model" . All operators are given on the level of dis- 
crete finite matrices. The eigenvalues of some "unper- 
turbed Hamiltonian" Hq form two "bands", i.e., a left 
and a right band, both of width Se, both with equidis- 
tant level spacing. Furthermore, there is a "perturba- 
tion" V consisting of transition operators, representing 
transitions between the two bands. The full Hamiltonian 
H = Ho + V reads 



Ho = ^^(56|*,i)(z,L| + V-^<5e|J,i?)(J 



j=0 



J=0 



(n — 1 \ ^ n — 1 

^ V,, KL)(j,i?|+H.c. ,A2 = - ^ It;. 



R\ 



2 



(3) 



where the Vij are chosen as random complex Gaussian 
numbers and A measures the perturbation strength. Note 
that the Vij are randomly chosen but fixed numbers in 
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this definition. This apphes for ah the numerics foUow- 
ing further below. \i,L) are the eigenstates of 

Hq which form the left (right) band. We only consider 
one relevant observable A having two eigenvalues, 1 and 
— 1, each rt-fold degenerate. Each eigenstate of A also 
corresponds to an eigenstate of i/o- The 1-eigenspace co- 
incides with the left band, the —1-eigenspace with the 
right band. If the expectation value of A is 1 (—1), the 
excitation probability is completely concentrated on the 
left (right) band and if the expectation value is 0, the ex- 
citation probability is equally partitioned between both 
bands. One may now define the projectors on the states 
forming the left (right) band as Pl (Pr)- Then A may 
also be expressed as A ^ Pl — Pr and therefore be in- 
terpreted as a "position observable" which measures the 
"occupation asymmetry" between the two bands. 



choice of the observable A, i.e., for other observables the 
dynamical behavior could be completely different. 
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FIG. 2: Dynamics of expectation value a{t). "Fluctuating" 

curve as generated from the Schrodinger Equation (dashed 

line) versus mean exponential behavior with the decay rate 

© R (solid line). System parameters: A*' = 6000, 5e = 0.75, 

V A = 0.00013. 



FIG. 1: Sketch of "quantum two site hopping model". The 
model consists of two energy bands, both of bandwidth Se, 
each containing n levels with equidistant level spacing. Tran- 
sitions between both bands are triggered by the "perturba- 
tion" V. 

In the following we analyze the dynamics of the ex- 
pectation value a{t) = {^lj{t)\A\^jj{t)) for pure states. 
Thereby we focus on situations, where the rough system 
parameters Se, n, A are chosen from a suitable parame- 
ter range, such that the dynamics of a{t), as generated 
by the full time-dependent Schrodinger Equation, result 
approximately as an exponential decay. The suitable pa- 
rameter range is essentially defined by 



IGtt^ n \^ 
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Eq. ^ implies that relaxation times have to be 
much longer than correlation times re, i.e., ^ tq 
with Tc « 47r/(5e, tr := 1/R and R « {tc cor- 

responds to the autocorrelation function of V.) Or in 
other words there has to be a seperation of time scales 
between tr and tc- So, if the criteria |4] are fulfilled, the 
expectation value dynamics are approximately described 
by a{t) a(0) exjp(—Rt) for the overwhelming majority 
of initial states [3 1 where the latter statement refers 
to the framework of "typicality" [T8l - l2l| . For a detailed 
discussion on the occurrence of exponential relaxation in 
this type of quantum models see An illustration 

of the graph of a{t) is given in Fig. [5] and may also be 
found in [13, [3 • The decay rate is consistent with a per- 
tinent projection operator approach or/and with Fermi's 
Golden Rule. 

Note that the exponential feature also depends on the 



However, the dynamics as obtained by exact diagonal- 
ization show deviations from an exact exponential decay 
("wiggles"), i.e., the complete dynamics may be viewed 
as a composition of a "regular" exponential part and 
some "irregular" looking part, the latter being perceived 
as fiuctuations of the expectation value a(t). 

3. Entropy Definition and Numerical Analysis of the 
FT As stated at the very beginning of there is a 
long and ongoing discussion about a microscopic non- 
equilibrium definition of entropy. Here we suggest an en- 
tropy S{a) for which eventually a transient FT as defined 
in Eq. ^ holds. S{a) is defined as S{a) = NS{p^cd{a)) 
[N — 2n being the dimension of the full Hilbert space). 
S{pi-cd{o-)) — — Tr{/3rcdln/9rod} corrcspouds to a "reduced" 
von Neumann-entropy with 



Prcd(a) = j^iil + a)PL 



(1 - a)PR) , 



(5) 



i.e., prcd is a "reduced" density matrix in the sense that 
it represents the maximum entropy state which is consis- 
tent with a given a. This of course is an implementation 
of Jaynes principle of maximum entropy as applied, e.g., 
in [231]. Eventually, S{a) exclusively depends on the ex- 
pectation value of our only observed quantity A. S{a) 
becomes maximal for a — 0, therefore a = is termed 
the "equilibrium value" in this context which is consis- 
tent with pertinent considerations within the framework 
of typicality 

The dynamics of a(t) now yield concrete entropy dy- 
namics, i.e., S(a(t)). We first numerically investigate 
whether these entropy dynamics fulfill a FT in the sense 
of Eq. ([T]). To this end we determine a large number 
of time-averaged entropy productions (I],- = {S{t + t) — 
S{t))/T) for different t and different trajectories, simply 
from numerically solving the Schrodinger equation. 
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The mean entropy production Eg is directly deter- 
mined from the decay rate R. Fig. [3] shows a distribu- 
tion of rescaled entropy productions, i.e., the mean value 
(as expected from the relaxation rate) and the variance 
(as expected by relation ([2])) are "scaled out", for a sys- 
tem with parameters TV = 6000, 5t = 0.75, A = 0.00013. 
Note that we use only one realization of the Hamilto- 
nian for these numerics, i.e., the "perturbation" ma- 
trix elements vij are randomly chosen but fixed num- 
bers. The data is taken from 1000 trajectories depart- 
ing from a set of initial states specified by a common 
a(0) but chosen at random otherwise. From each trajec- 
tory we calculate entropy productions over 10 segments 
which are located around an (average) expectation value 
of a « 0.2(0.21 - 0.16) at time t « 820 - 1050, i.e., there 
are 10000 entropy productions in total. The segment 
lengths r are chosen such that tc < r ^ r^j is fulfilled 
(i.e., r« 70, Tc« 15, Tj^ = 1177). 
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FIG. 3: Relative frequency distribution of normalized entropy 
productions ^ := (S — So)/-\/2Eo/r. The points represent 
the frequencies obtained from our numerics, the solid line is 
a Gaussian fit (/i « 0, cr « 0.94, cr^ ~ 0.88). System parame- 
ters: N = 6000, 5e = 0.75, A = 0.00013. 
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FIG. 4: Dependence of cr on a. a decreases with larger 
deviation from equilibrium (a = 0). System parameters: 
N = 6000, Se = 0.75, A = 0.00013. 



ponential relaxation of some quantum expectation value 
a(t) and an entropy definition as S{a) (see 3.), the rel- 
ative frequency of mean entropy productions over time 
steps which are larger than the correlation time but 
shorter than the relaxation time to fulfill a transient FT 
as given in Eq. (IT|). This is expected to hold within a 
regime around equilibrium the size of which may depend 
on the model. 

4- Stochastic Description of Expectation Value Dy- 
namics In the following we present some concepts of an 
(partially analytical) explanation for our main numeri- 
cal results. The key idea is to show that the dynamics 
of an expectation value a{t), as resulting from the fully 
deterministic Schrodinger Equation, may be described in 
terms of a time-discretized stochastic process of the fol- 
lowing Ornstein-Uhlenbeck-type: 



da, 



-Ra, 



^/2R/N■dw^ 



(6) 



Obviously the distribution of the entropy productions 
is well described by a Gaussian (see Gaussian fit in Fig. 
[3]). Wc find the mean and the width to be /i « and 
be cr = 0.94, respectively. The FT would have been ex- 
actly fulfilled if we had found ji — Q and a — 1. That 
is, the variance is slightly smaller than postulated by 
the FT. Recall that the analyzed trajectory segments 
are picked at times corresponding to expectation values 
of a w 0.2(0.21 — 0.16) which means significantly off- 
equilibrium (a = 0). Numerics (see Fig. 14]) suggest that 
the width gets closer to one with data taken closer to 
equilibrium for any model. The decrease of a for larger 
deviations from equilibrium cannot be fully explained 
within this framework. Further numerical investigations 
(which we omit here for brevity) indicate that the degree 
of the deviation from a — 1 when departing from equilib- 
rium seems to depend on the specific model. Note that 
the model at hand has not been designed to be realistic 
in any detail, it has only been designed to be simple. 

Based on these numerics we now formulate our first 
main result. We expect for models showing overall ex- 



Here the dwi are stochastic increments drawn from a 
standard Wiener process, i.e., the dwi are uncorrelated 
numbers drawn from a Gaussian distribution with zero 
mean, {dwi) — 0, and normalized variance, (dwidwj) = 
SijT. Again the description is meant to hold for r 
from a regime in between correlation and relaxation 
times. (Note that unlike in Nelson's approach [23| we 
do not argue here that a quantum probability density 
may be found from a Fokker-Planck equation for a cer- 
tain stochastic process, but that the evolution of a certain 
quantum expectation value has striking similarities with 
a sample path from some Ornstein-Uhlenbeck process.) 

By hypothesizing Eq. ([6]) and an entropy as given by 
a Taylor expansion of S{a) (see 3.) up to the quadratic 
order around the equilibrium value a = 0, i.e., S{a) sa 
A^ln(A^) — 1/2 Na^, one can analytically show that the 
probability distribution of entropy productions P(E,-) is 
Gaussian as well and fulfills the FT ([1]). Applying Ito's 
lemma the entropy, which is a function of the expectation 
value, also becomes a stochastic process that is described 



4 



by 



dS, = 2R{N\n{N)-S^ 



\)t- ^AR{N\^{N)- S^)-dw, 



(7) 

For large enough systems and large enough devia- 
tions from equilibrium, which is fulfilled in our setup, 
iVln(A^) - 5, is large compared to 1/2. (Both A^ln(iV) 
and Si become large for large N .) Consequently, the 
term 1/2 in the first expression on the r.h.s. of ([7]) may 
be neglected. Dividing by r and further introducing the 
abbreviations := dSi/r and Sq := 2R{N\n{N) - Si) 
one arrives at 



Ei — Eo 



(8) 



for a system of A'' = 6000, 5e = 0.75, A = 0. 00013 in di- 
cate that the distribution of the D[tp,t,T)^ N /2Rt as 
generated from evaluating 10000 segments along a sin- 
gle trajectory directly from the Schrodinger Equation, 
is approximately Gaussian with zero mean and ct^ « 1 
(Fig. [U (r « 70, Tc ~ 15, tt? — 1177, again we use only 
one realization of H, i.e., the Vij are fixed). This re- 
sult appears to be robust against variation of the model 
parameters and time step lengths t within the above 
regime, (iii) Moreover, Fig. [6] shows that the correla- 
tions < D[tp, t, t)D{'iP, t + At, r) >, as obtained from the 
same data that has been exploited to clarify the previ- 
ous issue (ii) , are nearly for the above chosen time step 
length T, such that the !)(■!/', i,T) can be considered as 
uncorrelated. 



where corresponds to the entropy production over the 
time step r and Eq is the mean entropy production. That 
is, the distribution of entropy productions is Gaussian 
and the mean value and the variance obey the relation 
given in which means that the FT ((IJ is fulfilled. 
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FIG. 5: Relative frequency distribution of normalized expec- 
tation value changes x •= D{ip,t,T)/ yj2RT/N . The points 
represent the frequencies obtained from our numerics, the 
solid line is a Gaussian fit (^ ~ 0, cr Ri 0.96, ~ 0.92). 
System parameters: = 6000, 5e = 0.75, A = 0.00013. 

We now analyze whether the dynamics of {ip\A{t)\4>) 
according to the Schrodinger Equation are consistent 
with the crucial features of ©. We especially regard 
the following points: (i) Is the overall behavior described 
by an exponential decay with the decay rate R7 (ii) Is 
the distribution of the quantity 

Di^, t, t) := (V|A(t + rm - (V|^(t)|^) (1 - Rr) (9) 

Gaussian with zero mean and variance ~ 2Rt/N1 
(Note that the deterministic D{ip,t,T) correspond to 
y/2R/N ■ dw^.) (in) Are D{'ip,t,T) and D{^,t + At,r) 
{At > t) uncorrelated? 

We start by investigating those points numerically for 
our model, (i) Here we simply resort to the results of 
[itI I , [3| ■ They give evidence that typically , up to small 
fiuctuations, the dynamics of {i)\A{t)\il}) is an exponen- 
tial decay with a rate R as given in 2. (ii) Numerics 
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FIG. 6: Correlations of K(At) =< D{tlj,t,T)D{i(),t + At,T) > 
/{2Rr /N). Each point is averaged over 1000 trajectory seg- 
ments, T « 70. System parameters: A'' — 6000, 5e — 0.75, 
A = 0.00013. 

Thus, the numerics may suggest that for < r <C 
Tji the deterministic D{ip,t,T) may be interpreted as 
stochastic variations yi and that a description by (|6]) is 
valid in and also close to equilibrium. 

In addition to the numerics we now provide some an- 
alytical considerations on the validity of a description as 
given by ([6]), which also address the points (i)-(iii). Note 
that these considerations are independent of any details 
of the model, except for exponential approach to equilib- 
rium. The statements are mathematically rigid, but on 
the level of ensemble arguments, e.g., it is shown below 
that the variance of the D{ip,t,T) is given by 2RtN^^ 
if the set (ensemble) of all possible states V' is consid- 
ered. Strictly speaking this holds in itself no implication 
on the distribution of the D{ip, t, r) encountered along a 
single trajectory, i.e., with respect to a set of t. Such a 
connection only arises from some sort of ergodicity. But 
following the ideas of von Neumann [2^, this ergodicity 
may be expected here. The literature provides a for- 
mula for the ensemble or "Hilbert space" average (HA) 
of the product of two quantum ex pectation values. The 
ensemble includes all possible pure states from a finite 
dimensional Hilbert space [26| : 

HA((V'|A|V')(^|B|V')) = {N{N + l))-iTr{AB} (10) 
for A,B traceless i.e., Tr{A} = Tr{B} = 0. Thus to 
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address (i) we simply compute 

HA((^|A(t + rm{ij\Am)) = ^^^^ , (11) 

Since we focus here on models featuring exponential de- 
cay we have Tr{A(T)A} cx e~^'^ « 1 — Rt. This justi- 
fies the usage of i? in (jS]). To address (ii) we calculate 
the Hilbert space average of D'^{il;,t,T) from ([9]), which 
should be a measure for the fluctuations. By means of 
Eqs. pUj) and (ITT|) one obtains 

IlA{D'ii;,t,T))) = (Tr{A2(t + r)} (12) 
- 2{1 - RT)Tr{A{t + T)A{t)} 
+ {l-RTfTr{AHt)}){N{N + l)y^ . 

Assuming the above exponential decay of the correlation 
function one obtains to linear order in r 

}iA{D^{^,t,T)) « 2RtN-^ , (13) 

i.e., the fluctuations become small for large N and are 
proportional to the relaxation rate. The result of (|13l) is 
exactly what we assumed as the width of the stochastic 
fluctuations in Eq. Regarding (iii), one finds by 

realizing that 

D{^, t, 2r) = t) + D{ijj, t + T,T) , (14) 

and exploiting (fT3)) that the HA of the correlations be- 
tween subsequent trajectory segments vanishes 

YiK{D{'iP,t,T)D{i!,t + T,T)) ^{) . (15) 



(ITSt may be straightforwardly extended to non- 
subsequent trajectory segments. 

To repeat, the above HA calculations address averages 
over all possible states in Hilbert space and therefore 
yield results only on the average behavior of all possible 
trajectories, but not directly for individual trajectories. 
Since Eq. ^ is certainly meant to describe individual 
trajectories, the HA results cannot strictly prove the cor- 
responding features of (jH). Nevertheless, the HA results 
are at least consistent with ([5]), which is not self-evident, 
and may in so far strengthen the trust in the validity of 
a description as given by 

5. Conclusion In this Letter we present an analysis 
of temporal fiuctuations in closed quantum systems fea- 
turing exponential relaxation of some expectation value. 
Numerics indicate that the dynamics of an entropy-like 
function are correctly described by a FT. This feature 
directly arises from Schrodinger-type dynamics, i.e., no 
assumption of, e.g., a Quantum Master Equation is in- 
volved. This result appears to be based on the finding 
that the underlying expectation value relaxation dynam- 
ics may be interpreted in terms of a stochastic process. 
In how far these results may be extended to quantities, 
which are better accessible by experiments, like, e.g., 
work, is an open question, particularly because work is 
not a quantum observable |27| . 

We thank C. Maes, A. Engel, H. Niemeyer and P. Maas 
for fruitful discussions. 
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